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Method of analyzing a data set comprising a tubular structure 



The present invention relates to a method of analyzing an object data set which 
comprises points in a multi-dimensional space and, in which a tubular structure occurs, said 
method comprising the following steps: 
a. choosing a starting position in or near the tubular structure; 
5 b. deriving a cutting plane through the tubular structure at the starting position, 

c. determining a number of points forming part of the surface of the tubular structure in the 
vicinity of the starting position; and 

d. calculating a gradient to the surface for each of said points. 

1 0 The international patent application EPOO/09505 (PHN 1 7.678) of the same 

applicant (not yet published) relates to a method of the type mentioned above. This method 
generally relates to the analysis of a tubular structure in a multi-dimensional space. 
According to the method described in the cited international patent application EPOO/09505 a 
self-adjusting probe is defined for analysis of the object data set. The self-adjusting probe 

1 5 comprises a sphere and a plane through the center of the sphere. The sphere should be 
positioned such that the tubular structure at least partially intersects the sphere. The plane 
should be oriented so as to be orthogonal to the tubular structure. When oriented correctly, 
the self-adjusting probe enables semi-automatic shape extraction of tube-like geometry. 

Such an object data set represents one or more properties of the object to be 

20 examined. The object data set notably relates to the density distribution in the object to be 
examined; in that case the data values are the local density values of (a part of) the object to 
be examined. The data values may also relate, for example, to the distribution of the 
temperature or the magnetization in the object. The multi-dimensional space is usually the 
three-dimensional space. The data values then relate to a volume distribution of the relevant 

25 property, for example, the density distribution in a volume of the object to be examined. The 
multi-dimensional space may also be two-dimensional. In that case the data values relate to a 
distribution of the relevant property in a plane through the object, for example the density 
distribution in a cross-section through the object. 
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The object data set can be acquired in a variety of ways. The object data set 
notably relates to a patient to be examined. Such an object data set can be acquired by means 
of various techniques such as 3D X-ray rotational angiography, computed tomography, 
magnetic resonance imaging or magnetic resonance angiography. 

The method described in the cited international patent application EPOO/09505 
is particularly suitable for analyzing the structure of a blood vessel. Several physical 
characteristics of the blood vessel can be determined, such as the diameter thereof, which are 
crucial for diagnosis and a safe treatment of, for example, a stenosis or an aneurysm. 

Currently the calculation of the above mentioned important step of orienting 
the cutting plane is based on a selection of surface points lying inside the sphere at a small 
distance from the plane. This is described in the article "Semi- Automatic Shape Extraction 
from Tube-like Geometry" by J. Bruijns, published in B. Girod, G. Greiner, H. Niemann, H.- 
P. Seidel (Eds.), Vision Modelling and Visualization 2000, Proceedings, November 22-24, 
2000, Saarbruecken. In the example described in the cited article the surface of a blood vessel 
is represented by means of triangles (a known method which will be explained later). 
However, the known selection of surface points is rather arbitrary and may unintentionally 
include points associated with neighboring vessels. This directly leads to inaccuracies in the 
determination of the physical characteristics of the blood vessel under analysis. 

The invention aims to provide a method of the type mentioned above which 
allows for more accurate calculation of the physical characteristics of the tubular stoicture 
under analysis. 

To this end, the method according to the invention also comprises the 
characterizing steps of: 

e) determining for each point a vector from the center of the tubular structure to said point; 

f) determining the angle between said vector and the gradient at said point; 

g) adding said point to a selection of points if said angle is equal to or smaller than a 
predetermined ceiling value; 

h) using said selection of points to calculate an orientation for the cutting plane such that the 
direction thereof is as parallel as possible to the longitudinal axis of the tubular structure 
at the starting position, and 

i) repeatmg the steps a) through h) for a new starting position along the tubular structure if 
necessary. 
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A preferred version of the method also comprises the steps of: defining a 
sphere which is at least partially intersected by the tubular structure, and performing the steps 
e) through g) only for surface points lying inside the sphere. When the radius of said sphere is 
chosen wisely, for example, slightly larger than the expected radius of the tubular structure 
under analysis, points within said sphere most likely belong to said tubular structure whereas 
points outside said sphere probably do not. Furthermore, the reduction of the number of 
points taken into account during the calculation also increases the calculation speed. 

In a further preferred version the steps e) through g) of the method are 
performed only for points lying at a predetermined maximum distance from the cutting plane. 
This allows for a further advantageous reduction of the number of points involved in the 
calculation thus further increasing the calculation speed. 

The invention also relates to a computer program for carrying out the method 
according to the invention. 

The invention will be described in detail hereinafter with reference to the 
attached drawing, in which: 

Figure 1 shows a cross section throu^ two adjacent tubular structures under 
analysis by way of the method according to the mvention. 

Figure 1 shows a cross section through two adjacent tubular structures that is, 
vessels 1 and 2. A sphere 3 is defined around the vessel 1. The center of the sphere 3 and the 
vessel 1 coincide and is indicated as c. The sphere 3 also partly intersects the vessel 2. 

The preferred version of the method according to the invention comprises the 
following steps. First a starting position is chosen in or near the tubular structure and a 
cutting plane is derived through the tubular structure at the starting position. These steps are 
also described in more detail in the international patent application EPOO/09505 which 
document is incorporated completely herein by way of reference. 

Next a number of points that form part of the surface of the tubular structure 
are defined in the vicinity of the starting position. In order to illustrate this step three points 
(also referred to as vertices) on the surfaces of the vessels are indicated. Point 4 is a surface 
point on the vessel 1 . Points 5 and 6 are surface points on the vessel 2. In the relevant art 
several methods are known to define the surface of a vessel in an object data set. When the 
object data set comprises a volumetric representation of the object by means of voxels, the 
surface of the vessel is defined by a set of boundary vessel voxels. This voxel set can be 
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determined by first separating the vessel voxels from the tissue voxels by means of knovm 
techniques such as a 'region growing algorithm'. Next it is determined which of the vessel 
voxels lie adjacent to tissue voxels. Such voxels are referred to as boundary vessel voxels. 
Preferably this determination is performed by finding the facing neighbors, that is neighbors 
5 which have one voxel face in common. As an alternative the surface of the vessel can be 
represented by surface triangles which can be generated by means of the so-called marching 
cubes algorithm which is known in the art. The marching cubes algorithm is described, for 
example, in the article: "Marching Cubes: A High Resolution 3D Surface Construction 
Algorithm", by Lorensen and Cline, Computer Graphics, Vol. 21, No 4, July 1987, 

10 In a next step gradient is calculated for each of the vessel surface points a 

gradient. This is performed by means of standard techniques which are among others 
described in the international patent application EPOO/09505, The resultant gradients for the 
points 4, 5 and 6 are indicated as v4, v5 and v6, respectively. 

In order to determine which of the three points should be included in the 

1 5 selection of points to be used for the calculation of the desired orthogonal orientation of the 
cutting plane, the method comprises the following characterizing steps. 

First for each point a vector is determined from the center of the tubular 
structure to said point. For the point 4 said vector is designated as vc4. For the point 5 s£iid 
vector is designated as vc5 and for the point 6 said vector is designated as vc6. 

20 Next for each surface point the angle is deteraiined between each of said 

vectors and the gradient at said point. If said angle is equal to or smaller than a predetermined 
ceiling value, said surface point is added to the above-mentioned selection of surface points 
to be used to orient the cutting plane. A ceiling value can be set by a user or can be 
automatically set. It has been found from experiments that a ceiling value of between 50 and 

25 70 degrees, for example approximately 60 degrees, gives reasonable results in the case of the 
current example relating to shape extraction for a blood vessel. 

Returning now to figure 1 some calculations will be performed in order to 
illustrate the method according to the invention. Starting with the point 4, the angle between 
v4 and vc4 is zero degrees. Assimie that the ceiling value is chosen to be 60 degrees. The 

30 point 4 is then added to the selection. This is rightly so, because the point 4 is a surface point 
on the vessel 1, that is, the vessel under analysis. However the angle between v5 and vc5 is 
more than 90 degrees. The point 5 will, therefore, be excluded from the selection under the 
current conditions. Again this is correct, since the point 5 is a surface point of the vessel 2, 
that is, the vessel neighboring the vessel investigated. 
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It is to be noted that the angle between v6 and vc6 appears to be smaller than 
60 degrees. The point 6 would unintentionally be added to the selection, since the point 6 is a 
surface point of the vessel 2. This can be prevented by adding another criterion relating to the 
selection of vertices, A useful second criterion is that only surface points lying inside the 
5 sphere 3 should be included in the selection. From figure 1 it is clear that the point 6 does not 
meet this criterion and would, therefore not be included in the selection after all. 

It is to be noted that a useful third criterion would be to perform the steps e) 
through g) for points lying at a predetermined maximum distance from the cutting plane. 
When de object data set comprises surface triangles, a ceiling value for this distance can be 
10 derived, for example, fi-om the length of the triangle edges (see also the above mentioned 
article by Bmijns). 

The order in which the three criteria of the method according to the invention 

are implemented may be varied. 

The selection of surface points thus created is then used to orient the cutting 
1 5 plane such that the direction thereof is as parallel as possible to the longitudinal axis of the 

tubular structure at the starting position. An example of this calculation is described in the 

article mentioned above, this is document is incorporated entirely herein by way of reference. 

The above steps of the method according to the invention can be repeated for each new 

starting position along the tubular structure if necessary. 
20 Now that the method of the invention has been explained, a skilled person will 

be able to translate the steps of the method into a computer program for carrying out the 

method. 

Summarizing the invention refers to a method for analyzing an object data set 
in a multi-dimensional space comprising a tubular structure. This method is especially useful 

25 in the field of medical diagnostics and treatment where the object notably is a patient to be 
examined. All kinds of tubular structures can be analyzed, such as blood vessels, but also 
brain cell or neuron cell structures. In analyzing blood vessels the method according to the 
invention allows for accurate determination of the size and diameter of a blood vessel which 
is essential for correct diagnosis of, for example, a stenosis or an aneurysm and for the safe 

3 0 treatment thereof. 

The invention is of course not limited to the version described or shown, but 
generally extends to any version that is within the scope of the appended claims as seen in the 
light of the foregoing description and drawings. 



